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Fernando Pessoa - 
“Sometimes I invent a metaphysics. With all the respectful scrupulousness of attention of someone 
engaged in real scientific work, it reaches the point where I may really be doing just that. The 
important thing is not to pride myself too much on all this, for pride is prejudicial to the precise 
impartiality of scientific exactitude." 


I'm not quite sure of the contents here yet. Nor when I will finish this. I never know how or where to begin, 
and I want the reader to know; that I have sat myself here to write without a clear idea of the language or 
events I will use to convey what it is that I want to express. And this meta-introduction, in all its immaturity 
speaks to what I will present—in both its immaturity and its honesty. Writing for myself is a fleeting thing, but 
I carry with me a very well-established and clear intention: to present an ontological and complete body of 
work that captures what I have learned throughout this. What I have learned and what I will depart to put 
forward are not answers; rather, my understanding of this, what this is. And in it, my own primal 
understanding nonetheless grows, and, in itself, forms a part of my current understanding. 


Both this, understanding, and their isomorphisms, I do now think to, in a deeply personal sense, ideas which 
may extend to a form of Universal or a form of Total. Whether by that I mean personally or collectively, 
whether what I may find to be universal may be merely relatable to whoever reads this; or something 
indigestible; or grotesque to another. And to that I am content. And what this is: I mean, and I do mean, my 
personal understanding, in my historical, epistemological, algebraic, and dialectical sense. 


I will attempt to present— what this writing would inherently lead to, if captured— Personally an ontological 
body of my phenomenology, a part of it; of my personal understanding; or conversely of a jumbled pretentious 
text. This all may likely stand incomplete, unfinished, and in both an acceptance of pretentiousness I feel a 
faith in spending my time towards, again; at least for these fleeting moments, time towards it all. 


I hate that word things. And maybe passion too, maybe I carry the idea of being passionate, when often I find 
myself too tired or stoned, high or bored, without means or lazy to, but perhaps it is that I have not constructed 
either the language or the actions, to perform such the task or thing that I am passionate about. Defining 
passion too, to a level that will satisfy what I am attempting to get at, I find difficult itself. 


And to who I am to write such a thing. The room I now live in is different to the one where I began writing 
this. I only write alone in my room. The room in which I sit is two and a half by four metres, and on the 
second floor of a contemporary-modern apartment: in the kind of bleak, corporate, real-estate, sense. But it's 
one that's horizontally clad with an acceptable false wood that mimics the Victorian New Zealand architecture 
in the surrounding suburbs. Stuck on the opposite ends of the symmetrical windows that run single file across 
the apartments of the second floor, the highest floor and the one in which I sit, are false French window 
shades. The ground level encircled by concrete stilts, to the roadside, the north-western side, a large pin-coded, 
metal-barred, roller gate and door. The ground floor occupied is parking spaces and an enclosed staircase on 
the left-hand side looking in from the street lead to 4 apartments above. The first floor of each apartment: an 
open space dining room and longue with a sliding door at the north-western end that leads out onto very small 
balcony, with wooden dividers that separate each apartment's balcony. Opposite, to the south-east, the door, 
and beside to the south-south-east the kitchenette stands tucked into the corner with a shiny, metallic newly 
bought modern appliances. Wooden surfaces lit by yellow LED lights, the temperature of which I'm unsure. A 
small a staircase somewhat hidden beside the kitchenette, slightly less deep into the room, leads to the two 
rooms and bathroom upstairs, my window at the southeastern end faces a courtyard with a building that 
mirrors mine across. The sides of the courtyard: two more apartment complexes, one of wood, whose rooms 
are tall and thin, with large windows, and a loft space I can see from my room on the right-hand side. To the 
north, my left-hand side, a concrete, three-storied, brick wall painted a blueish white. The apartment was 
rented unfurnished. I have recently bought: an Iranian hand carved wooden one-seater with sewn cushions that 
depict a girl and a deer, in what is quite an Art-Noveau style. Its arms: carved wooden lions. It is quite 
centrally in the dining room-lounge on the north-eastern side beside a wooden dining table and the table's 
accompanying two wooden chairs; A large metal coffee table that sits closest to the sliding balcony doors, 
itself with a glass top spray-painted black underneath; A glass desk for my room that measures 1.2m by 0.6m 
and sits on two metal A-framed legs in the south-south-east corner, it's legs two metal A's stuck together by 
two metal bars 0.6m across, and so two of each of those for each end of the table. With my personal computer 
tower, the machine I write this on, underneath the glass desk, equipped with an Intel i7 10700k, which my 
computer reads in this moment is running at 4.3 gigahertz, that is 430000000 cycles per second. 16GBs of 


random-access memory, an RTX 2070 GPU with a couple cores dedicated to matrix multiplication and 3TBs 
of a mix of solid state and disk space, of which 10085 remain. An “Old time", as the second hand furniture 
store put it, chair, with a curved 1080p 27” Samsung monitor on eye saver mode in front of me, a pair of KRK 
monitor speakers almost parallel but slight facing inward at the far corners of the table, an IBM Model M 
keyboard, Fernando Pessoa’s Book of Disquietude in front of me tucked under the monitor; beside a penguin 
collection of Allen Ginsberg’s Howl, Kaddish, and other poems. Tucked between the PC Monitor and the 
speakers are the books: Fyodor Dostoyevsky stories printed by Progress Publishers in Moscow, A stolen copy 
of the Feynman Lectures of Physics Vol: 3, Quantum Mechanics from the University of Auckland library 
(never to be returned), Theodor Adorno’s Minima Moralia, W B Yeats, and a photobook of myself gifted to 
me by my grandparents, Infinite Jest, half-way off off the Southern end of the table, the book cover of Martin 
Margiela’s women’s collections 1989-2009 as my mouse-pad for my RGB Hyper-X mouse, the book itself to 
my right side, with my weed grinder and phone above it, bong on the ground. An empty can of Coke Classic, 
Tom Ford Noir, an empty glass of kratom, and an empty wine glass, on the ground beside me an Uber Eats 
Bag, a Glow-in-the-dark Louis Vuitton styled large glass bong, a broken studded Versace belt, a New-Era cap 
above the speaker, and my bed on the northern side of the room. The staircase and second floor are carpeted a 
shallow but comfortable carpet, the dining room-lounge a fake wood, most likely a more chemically modern 
counterpart to linoleum. The walls, white. 


There are things that I do, and in moments, I find myself loving what I do. But why I do them: it's not that Iam 
unaware of the genetic, traumatic, and quite universally, probabilistically deterministic forces that position me 
into such roles. It's whether the reader wants to accept this as passion itself. 


Fernando Pessoa - 
"And naturally it is to find it necessary to reduce the spirit to some sort of physical matter surrounded 
by a space in which it can exist. All this depends on a great refinement of our inner feelings which, 
taken to its limit will no doubt reveal or create in us a genuine space like the space in which physical 
things exist but which in fact does not itself exist as a thing." 


I want to talk about structure, complexity and irreducible representations. It is to me unclear whether people’s 
writing about such things are isomorphic to what that I interpret those very things to mean. And so naturally as 
to whether I’m talking about mathematics yet, I already talk about the isometry of ideas, and here I’m talking 
about an equivalence of ideas, or at the very least, the relationships and mappings between them. 


Poincare - 
“How much a well-chosen word can economize thought. 


When our language is well chosen, it is astonishing to see how all the demonstrations made upon some 
known fact immediately become applicable to many new facts. Nothing has to be changed, not even 
the words, since the names are the same in the new cases. 


A word well-chosen very often causes the disappearance of exceptions to rules as announced in their 
former forms; it was for this purpose that the terms ‘negative quantities’, ‘imaginary quantities’, 
‘infinite points’, have been invented. And let us not forget that these exceptions are pernicious, for 
they conceal laws. Very well then, one of those marks by which we recognize the pregnancy of a 
result is in that it permits a happy innovation in our language. The mere fact is oftentimes without 
interest; it has been noted many times but has rendered no service to science; it becomes of value only 
on that day when some happily advised thinker perceives a relationship, which he indicates and 
symbolizes by a word. 


The physicists also do it just the same way. They invented the term ‘energy’, a word of very great 
fertility, because through the elimination of exceptions it established a law; because it gave the same 
name to things different in substance, but alike in form. 


How often mathematicians of the past must have run across groups without recognizing them and 
how, believing these groups to be isolated things, they have found them to be in close relationship 
without knowing why. Today we would say that they were looking right in the face of isomorphic 
groups. We feel now that in a group the substance interests us but very little; it is the form alone which 
matters, and so, when we once know well a single group, then we know through it all the isomorphic 


groups; thanks to the words ‘groups’ and ‘isomorphism’, which sum in a few syllables this subtle law 
and make it at once familiar to us all, we take our step at once and in so doing economize all effort of 
thought.” 


Algebra is the manipulation of symbols without (necessarily) regard for their meaning, especially in a way that 
can be formalized in logic, often cartesian logic. 


We talk about objects, about numbers, the differences between them, how they are unique, the relationships 
that they share with each other, how equal and unequal. When we are talking about the counting numbers, or 
natural numbers N = {1,2,3 ...} we are talking about the countably infinite set of positive whole numbers, each 
of the elements of this set have the same quality that they are all whole, positive, numbers where n + 1 always 
exists. The natural numbers form a totally ordered set under the standard "less than" relation. This means that 
for any two natural numbers a and b, one and only one of the following is true: a < b,a = b, ora > b. The 
natural numbers are a countably infinite totally ordered set. 


George Cantor – 
An aggregate (set) we are to understand any collection into a whole M of definite and separate objects 
of our intuition or our thought. The objects are called the elements of M. In signs we express x is an 
element of A x € A. An arbitrary set M made of components 1, 2 and 3 is given M = (1,2,3) We 
denote the uniting of many aggregates, M, N, P which have no common elements into a single 
aggregate (M,N,P), М = (1,2,3), N = (4,5,6),P = (7,8,9), Q = (M,N,P) =M UNUP= 
(1,2,3,4,5,6,7,8,9) . The elements of this aggregate are then the elements of M, of N, of P taken 
together. 


Bertrand Russel 
When mathematicians deal with what they call a manifold, aggregate, Menge, ensemble, or some 
equivalent name, it is common, especially where the number of terms involved is finite, to regard the 
object in question (which is in fact a class) as defined by the enumeration of its terms, and as 
consisting possibly of a single term, which in that case is the class. 


A Euler diagram of a set, 


The collection of all algebraic structures of a given type will usually be a proper class. Examples include the 
class of all groups, the class of all vector spaces, and many others. The paradoxes of naive set theory can be 
explained in terms of the inconsistent tacit assumption that "all classes are sets"., one could, for example, 
define a class of all classes that do not contain themselves, which would lead to a Russell paradox for classes. 


A partial order on a set is an arrangement such that certain pairs, for certain pairs of elements, one precedes the 
other, the word partial is used to indicate that not every pair of elements needs to be comparable. 


The elements of the power set of [x, у, z] ordered with respect to inclusion. 


Partial orders thus generalize total orders, in which every pair is comparable. A total order is some binary 
relation < on some set X, which satisfies the following for all a, b and c in X. 


a € a (reflexive) 
If a€bandb € cthena < c (transistive) 
Тј а < ђапађ <a епа = b (antisymmetric) 
Та < Ь от Ь <a (strongly connected) 


A well order on a set S is a total ordering on S with the property that every non-empty subset of S has a least 
element in this ordering. Let (P, €) be a preordered set and let S be a subset of P. That is S € Р. An element 
l € P is said to be the least element of S if l < 5 for all s € S. Every well-ordered set is uniquely order 
isomorphic to a unique ordinal number, called the order type of the well-ordered set. The well-ordering 
theorem, which is equivalent to the axiom of choice, states that every set can be well ordered. If a set is well 
ordered (or even if it merely admits a well-founded relation), the proof technique of transfinite induction can 
be used to prove that a given statement is true for all elements of the set. 


The integers Z = (..., 3, —2, —1,0,1,2,3, ... } can be formally constructed as the equivalence classes of 
ordered pairs of natural numbers (a, b), € N, with the equivalence relation ~ with the following rule 


(a, b)-(c, d) whena+d=b+c 


Addition and multiplication on integers can be defined in terms of equivalent operations on the natural 
numbers by using [(a, b)] to denote the equivalence class having (a, b) as a member one has: 


[(а, b)] + [Cc, d)] = [Ca + c, b + d)] 
[(a, b)] · [(c, d)] == [(ac + bd, ad + bc)] 


The rational numbers Q = (5, | a,b € Z}, e.g. 5, 8 may be built as the equivalence class of ordered pairs of 


integers. More precisely let (Z x (Z\{0})) be the set (m, n) of integers such n = 0. An equivalence relation is 
defined on the set by, 


(ті,п1)-(тҙ,") <> тіп; = тупа 
Addition and multiplication сап be defined by the following rules: 


(m4,74) + (Mz, n2) = (тіп; + nm; nng) 
(та, па) · (Mz, n2) = (тут, nnz) 


A metric space is a set together with a notion of distance between its elements. To see the utility of different 
notions of to see the utility of different notions of distance, consider the surface of the Earth as a set of points. 
We can measure the distance between two such points by the length of the shortest path along the surface, "as 
the crow flies"; this is particularly useful for shipping and aviation. We can also measure the straight-line 
distance between two points through the Earth's interior; this notion is, for example, natural in seismology, 
since it roughly corresponds to the length of time it takes for seismic waves to travel between those two points. 


The notion of distance encoded by the metric space axioms has relatively few requirements. This generality 
gives metric spaces a lot of flexibility. At the same time, the notion is strong enough to encode many intuitive 


facts about what distance means. This means that general results about metric spaces can be applied in many 
different contexts. 


Like many fundamental mathematical concepts, the metric on a metric space can be interpreted in many 
different ways. A particular metric may not be best thought of as measuring physical distance, but, instead, as 
the cost of changing from one state to another or the degree of difference between two objects. 


Formally a metric space is an ordered pair (M, а), where M is a set and d is a metric on M, i.e. a function. 
Satisfying the following axioms for all points x, y, z € M. The distance from a point to itself is zero, d(x, x) = 
0. The distance between two distinct points is always positive: If x = у, then d(x, у) > 0. The distance from 
x to y is always the same as the distance from y to x d(x, y) = d(y, x). The triangle inequality holds such 
that, d(x,z) € d(x, y) + d(y,z), this is a natural property of both physical and metaphorical notions of 
distance: you can arrive at z from x by taking a detour through y but this will not make your journey shorter 
than the direct path. 

d(p, q} = (ai — pi)? + (4 — рә)? 


q 


Ф 
а(р, q) => 
P2 Р m 


Фф = р 


Р qı 


A metric space M is called complete (or a Cauchy space) if every Cauchy sequence of points in M has a limit 
also in M. Intuitively, a space is complete if there are no “points missing” from it (inside or at the boundary). 
For instance, the set of rational numbers is not complete because e.g. V2 is “missing” from it, even though one 
can construct a Cauchy sequence of rational numbers that converges to it. There are sequences of rational 
numbers that converge at the limit to irrational numbers; these are Cauchy sequences having no limit in Q. 


e.g. the sequence defined 


: : 3 17 : " . "m 
Consists of rational numbers (1, тез ); however, it converges to the irrational V2 which does not exist in 


0. 


It is always possible to “fill all the holes”, leading to the completion of a given space. 


3363 


1393 | 


239 985 | 2378 577 


= 169 | 408 


A Cauchy sequence is a sequence whose elements become arbitrarily close to each other as the sequence 
progresses. 


We will consider the set of all Cauchy sequences that are defined as follows: 
C = о, лем | Xn € Ф and (Xn) nen is a Cauchy sequence} 


For two elements, that is, for two sequences (а. ) nen, (bn) nen, define: 


(Qn) nen~ (hn) nen ‘= (an A bn) nen 


It is a convergent sequence with limit 0, ~ is an equivalence relation on C (reflexive, symmetric, transitive), an 
equivalence class exists such that 


(ха лемі = (Cas nen | (Xn)neN € Cj 


We can define the real numbers 


C 
R = — = (бољем) Сален € C) 
This is the set of all points on the number line. 


[Can)nen]~ + [On)nenl~ : [Can + bn)nenl~ (well — defined) 
(ап) лем] ~ [bn)nenl~ += [Can : bn) nen]~ (well — defined) 


[((Qn)nenle < [banen]. “> 36250 INEN Vnz2N:ó-c«b,-—a, 


A sequence X1, X2, X5, ... In a metric space (X, d) is called Cauchy if for every real number r > 0 there is a 
positive integer N such that for all positive integers m,n > N, d(Xm: Xn) € т. A metric space (X, d) is 
complete if any of the following equivalent conditions are satisfied. That every Cauchy sequence of points in 
X has a limit that is also in X. Every Cauchy sequence in X converges in X (that is, to some point in X). A 
distance metric for the real numbers is defined d(x, y) = |y — x| by the absolute difference. The absolute 
value is the non-negative value regardless of sign. 


Natu ral, N Start with the counting numbers (zero may be included). 
e. е е LÀ » 
0 1 2 3 
Integer, Z Extend the line backward to include the negatives. 
4 9 e. e. e e e e ь 
-3 -2 -і 0 1 2 3 
Rational, Q Insert all the fractions. 
-23/4 -114  -15 1/ 11/3 23/4 
4• + 6: • · 0 > . · о > о @Q@eee@eee о. .· .· о. op 
-3 -2 -1 0 1 2 3 
Real Algebraic, A, Insert all the roots. 
-У5 -v2 -1/; 1/ v2 v5 
C EELE: ELELE, ELELE: 0-0-0 -@- ELELE! ELELE? ELI рь 
-3 -2 -1 0 1 2 3 
Real, R Fill in all the numbers to make a continuous line. 
-1 -e -У2 -1/2 1/2 v2 е т 
CERERE SO CETERE Pee ee CEEE CEEE Pee “ . . Pee >» 
-3 -2 -1 0 1 2 3 


The cardinality of a set is a measure of the number of elements in a set, two sets A and B have the same 
cardinality if there exists a bijection (a.k.a., one-to-one correspondence) from A to B. Any set X with 
cardinality less than that of the natural numbers or |X| < |№ is said to be a finite set, e.g. a finite set X = 
{1,2,3} with 3 elements. Any set that has the same cardinality as the set of natural numbers, or |X| = |N| = 
No, is said to be countably infinite. Any set with cardinality greater than that of the natural numbers or |X| > 
INI is said to be uncountable. The real numbers R are uncountable, any continuous range in the real numbers 
is also infinite. e.g. the subset of real numbers X = {x | x € R, a € x < b} generated from the range (a, b) 
is infinite, whereas any subset of the natural numbers is finite. 


Transcen- 


dental 


Fernando Pessoa - 
"The gloriously infinitesimal remains where it is, never goes beyond what it is, and it lives free and 
independent. The useless and futile create intervals of humble static in our real lives. The mere 
insignificant existence of a pin stuck in a piece of ribbon provokes in my soul all manner of dreams 
and wondrous delights. I pity those who do not recognise the importance of such things." 


The natural conclusion that Caeiro draws from the smallness of his hometown is the one that I do, reading 
those words, receiving from them what I feel to be both inspiration and liberation. Because his village is small 
you can see more of the world there, than you can in the city, and in this sense the village is larger than the 
city. I am equal in size to whatever I see, not hemmed in by the size that I am. I am filled with the same 
conviction Pessoa was filled with, of their ability to recognise the heavens into new constellations, I am equal 
in size to whatever I see. 


I think of my present, or the present right in front of me as a point, or an infinitesimally small point in a space, 
with rigour a Hilbert space. At this point, it is the summation, it is the everything that everything bears down 
upon me at this point, I can think of the geodesics of light through time and space piercing through this point, 
raging from the sun, my computer screen. I will embark latter to describe the complexity of the very structure 
of what I am now loosely calling a point. In this infinitesimal point of space lies the most tangled complexity, 
the same complexity that knits itself across everything. 


Now is right where I am, and where I am, is everything now. Personally, I imagine two twin doves as I 
describe this, flying beside each other their whole life, eating the same, desiring the same, with the same bird 
song and same view each day. I point to physical separation, distance as the thing that separates our universal 
experience. And not just the distance between the pen and paper on my desk in a static moment in space, but 
the distance between the pen and paper on my desk across time. 


As to my description of space as a Hilbert space rather than a Euclidean space, just space, or another common 
description of three-dimensional space. It is for the rigour required for this such topic, a space in which things 
exist has a form in itself, and to begin to describe the mechanisms of space and time it is to create abstract 
spaces that have properties which are useful in their ability to help with the description of the former. 


Euclidean space is intended to represent physical space. In the description of a point in three-dimensional 
space, its coordinate: Cartesian and Polar systems are often used. Generally, a coordinate in a Euclidean space 
of any dimension, that is to be able to give a description of the location of a point on a line, or a point on a 
sheet of paper, or a point in the general three-dimensional space around us, the Cartesian coordinates are the 
signed— meaning being either a positive or negative number—, distance from a central point, along 
perpendicular axes. For a point on a line, that would be a single number, for that on a sheet of paper, formally 
a plane, two coordinates (x, y), where y is perpendicular to x. and so forth for higher dimensions, (x, y, z) 
where z is perpendicular to y. and etc. A general way to think of any Euclidean space is a set of points 
satisfying certain relationships, expressible in terms of distance and angle. 


Euclidean two-dimensional vector space V is equivalent to IR? = R x IR. The standard basis vectors of IR? are 
Х = ил = (1,0) and у = v, = (0,1). Any other number in the vector space сап be reached as some 
linear combination of the basis vectors. Any other pair of linear vectors in IR? also forms a basis in R? 


The Euclidean plane IR? can be equipped with many different metrics, The Euclidean distance ‚ог Euclidean 
metric is defined d; (Ол, y1), (x2, y2)) = V (x2 — х1)? + (у; — y1)? 


A point in 3-dimensional (Euclidean) vector space given in cartesian coordinates by v = (5,9,2) Є R? = 
RXRXR v= aî+bĵ+ ck with standard basis elements = (1,0,0), 9 = (0,1,0), 2 = (0,0,1), these 
basis elements are sometimes denoted i,j, k or eo, €1, 62. 


The set of real numbers also carries naturally the structure of a topological space and as such R is called the 
real number line. The notion of topological space aims to axiomatize the idea of a space as a collection of 
points that hang together (“cohere”) in a continuous way. Roughly speaking, a topology on a set “of points" 
prescribes which subsets are to be considered “neighbourhoods” of the points they contain. Various conditions 
or axioms must be satisfied in order for such neighbourhood systems to form a topology, but one of the most 
important is that for any two neighbourhoods of a point, their intersection must also be a neighbourhood of 
that point. Together with the cartesian product X - the cartesian spaces R”, for natural numbers n € N, form 
the classical spaces of physics in n dimensions. 


A manifold is a space that looks locally like a cartesian space, commonly a finite dimensional cartesian space 
R”, in which case one speaks of a manifold of dimension n or an n-fold, but possibly an infinite dimensional 
topological vector space, in which case is an infinite dimension manifold. A differentiable manifold is a 
topological space that is homeomorphic to a Euclidean space of some degree of differentiability. If one speaks 


of arbitrary or infinite differentiability, then one speaks of smooth manifolds. In any case, the type of geometry 
embodied in a particular flavour of manifold is controlled by a particular groupoid or, more generally, 
category of transformations which preserves whatever geometric features one is interested in. 


That is a topological space X is called locally Euclidean if there is a non-negative integer n such that every 
point in X has a neighbourhood which is homeomorphic to real n-space IR". The real coordinate space R” is 
an n-manifold, Any discrete space is a 0-dimensional manifold, A circle is a compact 1-manifold. Any space 
where there can be two distinct point x and y such that the neighbourhoods of x and y are unique is said to be 
a Hausdorff space. Almost all spaces encountered in analysis are Hausdorff; most importantly, the real 
numbers (under the standard metric topology on real numbers) are a Hausdorff space. More generally, all 
metric spaces are Hausdorff. 


A function from a set X to Y assigns each element of X exactly one element of Y. The set X is called the 
domain of the function, and the set Y is called the codomain. 


f:X Y 


e.g. [ Ко R, 
f(x) =x? +2 


A surjective is a function f such that, for every element y of the function's codomain, there exists at least one 
element x in the function's domain such that f(x) = y. In other words, for a function f : X — У, the codomain У 
is the image of the function's domain X, It is not required that x be unique; the function f may map one or 
more elements of X to the same element of Y. 


An injective function is a function f that maps distinct elements of its domain to distinct elements; that is, x1 > 
x2 implies f(x1) Z f(x2). (Equivalently, f(x1) = f(x2) implies x1 = x2 in the equivalent contrapositive 
statement.) In other words, every element of the function's codomain is the image of at most one element of its 
domain. The term one-to-one function must not be confused with one-to-one correspondence that refers to 
bijective functions, which are functions such that each element in the codomain is an image of exactly one 
element in the domain. 


A bijective function, or one-to-one correspondence between two mathematical sets is a function such that each 
element of the second set (the codomain) is mapped to from exactly one element of the first set (the domain). 
Equivalently, a bijection is a relation between two sets such that each element of either set is paired with 
exactly one element of the other set. 


A Homomorphism is an isomorphism if and only if it is bijective. In most cases bijective mappings are 
isomorphisms. 


As for an isomorphism, it is the same as the inverse morphism. The multiplicative identity of X operating on 
X gives X. And note I% -X = X , for example 1 - 2 = 2. Where 1 is the multiplicative identity of the real 
numbers. 

In all cases and using the composition of functions it is that f^! o f = I, , more simply we can say if a 
function operates on x and gives some number, then taking the inverse of that function and applying it on the 
number produced will give the original x. 


f (feo) = x, for every x eX 


7 
<1 


f 
P umm 
f 
в 
f! 


A morphism is a map between objects. Two objects x and y are isomorphic if there exists an isomorphism 
from x to y (or equivalently, from y to x). A morphism that has an inverse morphism is called an 
isomorphism. It is an abstraction of homomorphism, a function between two algebras that preserves the 
structure. This itself is even an abstraction of a more descript function that has some mapping from x to y i.e. 
f:x ә y. That is that x and y share a similar structure. 


x—.y 
g 


Z 


The function f is a morphism between X and Y, as g is from Y to 2, the composition of functions 2 and f gives a 
morphism from X to Z 


gof 


Given arbitrary functions, f (x) = x + 3 апа g(y) = у + 3, and an arbitrary set X = (1,2,3), we have, Y=f()=X+3 = 
(4,5,6), 2 = g(Y) = (45,6) + 3 = (7,8,9), 
Z = g 0 f =g(f(x)) = (025 + 3) + 3 = ((Х + 3) +3) = (7,8,9) 


In category theory, the idea of functions can be abstracted further into the idea of functors, a function is an 
element-wise mapping from X to Y , but it does not inherently preserve or map the structure of X onto Y. A 
set X, barren of additional structure is comprised solely of elements x € X. A function from one set to another 
merely sends elements from the first to the second. A category C is comprised of objects ob(C) and between 
any two objects A, B, a class of morphisms, denoted hom (A, B). The objects А, В need not be distinct; there is 
always the identity morphism defined id, € hom (А, A) so that f о id, = f and id4 ° g = g for any f € 
hom(A,-) and any g € һот(:, A). 


So, we can say that if f: X > Y is a homomorphism, then the function acting on the identity of X is equal to 
the identity of Y, and that the function acting on the identity of Y is equal to the identity of X. 


f Ux) = ly 
Proof : take, x ЄХ, 
f (xly) = Ғооғау) 
Қау) = ҒО) 


(FIFO) FG) = (EF 621 * £C) 
lyf (ly) =I, 
f (ly) = 1, 


Such that, 
X=Y 


There is a deep form of perseveration of structure, and to a degree what information can be held within that 
structure by the bijective mapping, the isomorphism. 
In some cases, 

X=Y,orX =Y 


We see that the groups are the same in a structural way, but their elements are different. Their elements may 
be equal. Equality is when two objects are exactly the same, and everything that is true about one object is true 
about the other, while an isomorphism implies everything that is true about a designated part of one object's 
structure is true about the other's 


In cases, in lower categories, we may assume an isomorphism to be a form of equivalence or equality between 
two objects, groups, numbers, algebras. 


And for this inverse morphism to exist, there must be underlying structure between X and Y such that they are 
the same. There is a need for good definitions. 


Bryan Garner - 
“What’s wrong with genteelism’s like ‘prior to’ and ‘subsequent to 

David Foster Wallace - 
“Well, I have trouble parsing your question. Genteelism seems to me to be an overly charitable way to 
characterize them. They are to me like puff words. They are like using “utilize” instead of “use”, 
which, in 99 cases out of 100, is just stupid. Or ‘individual’ for *person'. More syllables it's just 
puffed up. Why say ‘prior to’ rather than ‘before’? Everybody knows what ‘before’ means, it's fewer 
words. And I think, technically, given the Latin roots, it should be ‘posterior to’, if you're going to use 
‘prior to’. So, if you're saying ‘prior to’ and ‘subsequent to’, you are in fact, in a very high-level way, 
messing up grammatically. But would you ever want to say, ‘posterior to’? No. But you will notice, 
this is the downside of starting to pay attention, you start noticing , you know, all the people who say, 
‘at this time’, rather than, ‘now’. Why did they just take up one-third of a second of my lifetime 
making me parse, ‘at this time’, rather than just saying ‘now’, to me. And you start getting bugged. 
But you get to become more careful and attentive in your own writing, so you become an agent of 
light and goodness, rather than the evil that's all around." 
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Consider the following two strings of 32 lowercase letters and digits: 


Abababababababababababababababab, and 
4c1j5b2p0cv4w1x8rx2y39umgw5q85s7 


The first string has a short English-language description, namely "write ab 16 times", which consists of 17 
characters. The second one has no obvious simple description (using the same character set) other than 
writing down the string itself, i.e., "write 4clj5b2pOcv4w1x8rx2y39umgw5q85s7" which has 38 characters. 
Hence the operation of writing the first string can be said to have "less complexity" than writing the second. 


More formally, the complexity of a string is the length of the shortest possible description of the string in some 
fixed universal description language (the sensitivity of complexity relative to the choice of description 
language is discussed below). It can be shown that the Kolmogorov complexity of any string cannot be more 
than a few bytes larger than the length of the string itself. Strings like the "abab" example above, whose 
Kolmogorov complexity is small relative to the string's size, are not considered to be complex. The 
Kolmogorov complexity can be defined for any describable object. 


An irreducible representation is a representation that has no smaller non-trivial representations "sitting inside 
it”. Given some algebraic structure, such as a group, equipped with a notion of (linear) representation, an 
irreducible representation is a representation that has no nontrivial proper sub-object in the category of all 
representations in question and yet which is not itself trivial either. In other words, an irrep is a simple object 
in the category of representations. A representation that has proper nontrivial sub-representations but cannot be 


decomposed into a direct sum of such representations is an indecomposable representation but still reducible. 
In good cases for finite dimensional representations, the two notions (irreducible, indecomposable) coincide. 


Be the north wind, it's the one that hits here, the most, my home, and for the poetry of it ГЇЇ say, “where the 
virgin laid". For it's the poetry of this narrative, and that one that can be written about the algorithm, my life. I 
feel, as the nights were colder before. They arrive, cats, drawn by the yoke of bulls, and before they go, they 
offer everything of their gentle protection. For I know it's the cradle, of what is coming. The cats always come 
past, the past I was shown, on the steps where played with the black cat, in a place near a place named Nelson, 
New Zealand. I built things a lot when I was younger—rules for games, role-playing as a child—in relation to 
the things that I was, and am still, passionate about: Imaginatively in my backyard I would advance from the 
neolithic to the contemporary. On my personal computer I would do the same, progressing from the raw earth 
underneath the Robinia tree, through logic and the creation of language, into thoughts themselves, advancing 
through feudalism, castles, to imperialism and back to my cotemporary living room. When I entered 
preparatory school, rockets from PVC tubes, cooking stump remover (potassium nitrate) as an oxidizer, 
sometimes mixed with rust power, with white sugar, sucrose, as the fuel. The cooking process over a camp 
burner to form a solid candy like fuel instead of powder to be stored inside the rocket. I note this cooking 
process led to explosion twice. And it was either the metal end of a garden hose—the same New Zealand hose 
that would become smaller and smaller as the years went by—and once moulded concrete... which led to 
explosion upon ignition. Ignition occurred by connecting a 9V battery to a switch from my local circuitry store 
and nichrome out of a smashed filament bulb which would ignite outside of its vacuum once connected to a 
power source. This hobby evolved into turbine jet engine building in high school. The PVC tubes were 
replaced with stainless steel, and the turbines were cut from sheets that sent the blade of my small circular saw 
flying, returning me to the computer. This led me to pursue a science degree and eventually embark on post- 
graduate study in quantum physics. My passion for linear algebra in my undergraduate studies led me to this 
path, though I ultimately dropped out—a decision influenced by the logic I constructed from the language in 
my head." And ultimately in the pursuit of my passion, which I hate to declare out of the cringe declaration of 
practice, the pursuit to express my feelings and perception by whatever means. 


The modern concept of limit originates from Proposition X.1 of Euclid's Elements, which forms the basis of 
the Method of exhaustion found in Euclid and Archimedes: "Two unequal magnitudes being set out, if from 
the greater there is subtracted a magnitude greater than its half, and from that which is left a magnitude greater 
than its half, and if this process is repeated continually, then there will be left some magnitude less than the 
lesser magnitude set out." 


The expression “0.999...” should be interpreted as the limit of the sequence 0.9, 0.99, 0.999, ... and so on. 
This sequence can be rigorously shown to have the limit 1. And therefore, this expression is meaningfully 
interpreted as having the value of 1. Formally suppose ад, аз, ... is a sequence of real numbers. When the limit 
of the sequence exists, the real number L is the limit of the sequence if an only if for every real number є > 0, 
there exists a natural number N such that for all n > М we have |а, — L| < =. The common notation 


lim (a4) = L 
п- оо 
is read “the limit of а, as n approaches infinity equals L”. 


A differentiable function of one real variable function is a function whose derivative exists at each point in its 
domain. 


A function f : U > R, defined on the open set U C К is said to be differentiable at a € U if the derivative 


vs cs (Ғазю- ја) 
түше. 


exists. This implies the function is continuous at a. 


The derivative of a function of a single variable at a chosen input value, when it exists, is the slope of the 
tangent line to the graph of the function at that point. The tangent line is the best linear approximation of the 
function near that input value. For this reason, the derivative is often described as the instantaneous rate of 
change, the ratio of the instantaneous change in the dependent variable to that of the independent variable. The 
process of finding a derivative is called differentiation. 


Derivative of f is the rate of change of f 
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Euclid - 
"The curved line a is the first species of quantity, which has only one dimension, namely length, 
without any width nor depth, and is nothing else than the flow or run of the point which will leave 
from its imaginary moving some vestige in length, exempt of any width." 


This definition of a curve has been formalized in modern mathematics as: A curve is the image of an interval 


to a topological space by a continuous function. 
Co Neve хло S у солк КА 


Suppose а curve is given by some function f(x), a curve by definition is continuous such that it has а 
derivative. Suppose the path of some object across time is given by f(t). The derivative of f with respect to 
time is the velocity f'(t) = v(t), the derivative of the velocity gives acceleration f” (t) = v'(t) = a(t). Then 
that object path through time with respect to some origin is given by, 


БӨ = х(0) + y(Oj + 2(0) 
The instantaneous velocity is given by 
ај dx(t), dyt), dz(t) Ё 


QUT dt ^ dt^" dt 


The acceleration 
df Я) dv Е d?x(t) Й d? y(t) А d?z(t) ў 


г) Еа ааа -- рар салад 
a) dt? dt dt? ' dt? ) dt? 


1 w.r.t some origin, 2 with intrinsic coordinates 


For functions with multiple parameters say such as a function that changes independently in the x, y and z 
directions the derivative in respect to one of those directions is given by the partial derivative in that direction. 


e.g. for a two-dimensional function, i.e. a surface. f : x,y > R3, the derivative at a point (a,b) in the x 
direction is given by 
Of (a,b) 


Ox 


ð 
Line has slope 2 (a, b) 


Graph of f(x,b) 
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When we take the derivative of such a function we are taking the sum of partial derivatives, 


dF д д д 
dX Ox ду 02 

The partial derivative, that is the gradient or slope of a function in some defined direction, or variables is 

defined as a limit. Let U be an open subset of R” and f: U > Ra function. The partial derivative of f at the 


point a = (a4,a5,...,a4) € U with respect to the ith variable x; is defined as 


flat ће;) – f(a) 

h 
Where e; is the unit vector of the ith variable x;. 
An important example of a function of several variables is the case of a scalar valued function f (x1,..., Xn) on 
a domain in Euclidean space R” (e.g. on R2, f(x, у, z)) , in this case f has a partial derivative with respect to 
each variable (x, у, 2) or (x1,x2, x3) . At the point a, these partial derivatives define the vector 


a | 
qu TIN 


n : " 
Vf(a) = Gr (a),... „21. (а)), more generally we can abuse this notation and define, V = 5] ej — [-] e + 
n j= 


2.” a 
[бу 162 +. = Hazle 


i. The Circle 


Imagine any right-angle triangle, with its point, the one at the peak of its adjacent and hypotenuse, dead centre 
in a circle, and that the point at the peak formed by its opposite and hypotenuse lies on the circle. 

The angle that is formed between the opposite and hypotenuse, for any of these triangles we have described, is 
sind. 


We can describe any point on the circle then algebraically. 
x = т соѕ(0), у = т 5110, where x lies on some axis, and y lies on an axis, perpendicular, orthogonal, at 0 on 
the x axis, your standard graph. 


Length, magnitude, and Euclidean norm is given by | х I= 4 | x?, the length of a two dimensional vector or 
"arrow" is given by the square root of the sum of the square of its two components, by the Pythagorean 


theorem || x ||= J x? + y? in three dimensions || x ||= x? + y? + 22 etc... . Unitary matrices and unit vectors, 
unit complex numbers etc. all have length / norm of 1, unit length. 


|x| is the notation for the absolute operator forcing x value of x to be positive, and real. The norm is also 
always positive and real. 


51, is the 1-sphere and it is the circle, the circumference of a disk, in the two-dimensional plane. 
Homeomorphic to the real projective line. 
In two dimensional cartesian coordinates an element on the 1-sphere, i.e a point on the circumference of a 
circle with radius 1 can be parameterized by two coordinates, (x, y) such that ./x? + y? = 

i. The Circle. A New Approach. 


The imaginary unit i is defined solely by the property that its square is —1: 


With i defined this way it follows that i and —1 are both square roots of —1 


Given 
i? = —1 or Vi = У—1 
їх[=—1 
—1хї=—[ 
—ixi-1 
1xi-i 
More so 
i = 1, i! i, i = –1, i = -i, 
й =1, 9-4 ®=-1, %4--% 


That is the powers of i are cyclic. 
A complex number is an expression of the form a + ib 
The set of all complex numbers: 


C= {a +ib|a,bER, i? = –1) 


Im 


z=atbi 


The complex numbers C can be written as ordered pairs (a,b) of real numbers a and b, with the addition 
operator being component wise and with multiplication defined by 


(a, b)(c, d) = (ac — bd, ad + bc) 


The complex number whose second component is zero is associated with a real number: the complex number 
(a, 0) is associated with the real number а. 


The complex conjugate (a, b)* of (a,b) is given by (a,b)* = (a*, —b) = (a,—b), since a is a real number 
and is its own conjugate. The conjugate has the property that 


(a, b) (a, b) = (aa + bb, ab — ba) = (a? + b?,0) 


Which is a non-negative real number. In this way conjugate defines a norm, making the complex numbers a 


1 
normed vector space over the real numbers: the norm of a complex number 2 is |z| = (27272. 


Given a cartesian coordinate system, the unit n-sphere can be defined as 5" = {x € R"*?: || x ||= 1).where n is 
the dimensionality of the sphere. this reads as the n-sphere is equal to the set of all objects іп ІК"! with a 
magnitude equal to one. 


This is a very common structure in 2 and 3 dimensions, the circle and the sphere, the n-sphere generalizes the 
geometric structure of the sphere and extends it to any n dimensions. 


, it can be also thought of as a vector with magnitude, length of 1 parametrized by an angle 0, between 0 and 
360 degrees or ф € (0,27) radians around the origin, this is isomorphic to a unit complex number, given 
Euler's identity e" + 1 = 0 and e? = cos + ising . In the complex plane, which is a special interpretation of 
the common cartesian x-y plane, i is the point located one unit away from the origin along the imaginary axis 
which is orthogonal to the real axis. 


i.e. the circle group 5? is equivalent to the unit complex numbers Т, or the unit real numbers іп R?. 


81 = {x € R°: |I x 1= 1} = T = {z € C: |I z I= 1}, 


e? = cos o + i sin ф 


ry 


As a complex number consists of two independent real numbers, they form a two-dimensional vector space 
over the real numbers. Besides being of higher dimension, the complex numbers can be said to lack one 
algebraic property of the real numbers. Every Dedekind complete totally ordered field is isomorphic to the 
reals, where squares are necessarily non-negative in an ordered field. This implies that the complex numbers 
cannot be ordered since the square of the imaginary unit i is —1. 


An example of a geometric series found by Archimedes. 
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Wherea = “andr = 14, 


1 
ae 
That is, 
1 со 
=) me 1+х+х?+х% Gee 
1-х 
n=0 
which is valid for all |x| < 1 
A Maclaurin series is an infinite series of the form: 
Ginx” = dg + ах + ах? + 
n=0 
3 5 7 
__ 7 2 T cos(x) = 1— — + — – — + 
пл. =a a 
со TL 29 
Е У) (-1) NE = У. 
n=0 (2n + 1)! n=0 


And we see the similarity in the taylor expansion to that of sine and cosine, in fact we see their equality: 


2=а+ір 
со 
E 7 7° _ 2" 
е“ = titat = Рт] 


п=0 
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= cosx + isin 2, 


The group of rotations in 2D is SO(2) = S+, that is the orthogonal group O(n) being the group of distance 
preserving transformations of an n-dimensional space and the special orthogonal group SO(n) which contains 


the identity element. SO(2) is isomorphic to the group U(1) of complex numbers of unit length. 


Multiplication on the circle group is equivalent to the addition of angles, therefore it is the multiplicative group 
of all complex numbers with an absolute value (consider this intuitively as length or radius) of 1, simply the 
unit complex numbers. 


cos@ —sin@ 


is the group isomorphism. 
sin@ 17 group р 


This follows immediately from the fact the mapping e/? ^ ( 
Geometrically observe that U(1) is the unit circle St (50(2) = U(1) = 51). If you view the real and complex 
numbers purely as sets ignoring all possible structure on them, then the elements of R? and C are isomorphic in 
the category of sets due to the bijection (a, b) e а + bi. In real vector spaces R x R has the basis [(1,0), (0,1)] 
or (x, y}and C has a (real) basis (1, i}, since these real vector spaces both have dimension 2, they are 
isomorphic. Every plane rotation p, by an angle 0 is represented by multiplication by the complex number 

e? € U(1), in that sense for any two complex numbers 2,2' € С 

герб z = ez. 


Integration 


In the one dimensional case, with respect to integration in the simplest terms we want to find the area under 
some curve, from an interval (a through to b), sometimes we want to find the area under a curve that interval 
is infinite, and sometimes we want to find the area under a curve whose parts are both positive and negative, 
for all that means, “negative area". It is also thought of as an antiderivative, Velocity is the integral to 
acceleration, it’s the antiderivative to acceleration, for if acceleration represents the rate at which velocity is 
changing at each point in time, then its velocity is the sum of its acceleration across all points in time, or 
between some interval. 


“I have to pay a certain sum, which I have collected in my pocket. I take the bills and coins out of my pocket 
and give them to the creditor in the order I find them until I have reached the total sum. This is the Riemann 
integral. But I can proceed differently. After I have taken all the money out of my pocket, I order the bills and 
coins according to identical values and then I pay the several heaps one after the other to the creditor. This is 
my integral.” (letter from Henri Lebesgue to Paul Montel) 


The Riemann sum, let f be a real-valued function defined on the interval (а, b] The Riemann sum of f with 
respect to the tagged partition xo, ..., Ху together with to,...,t,-1 is 


n-1 
У б x) 
i=0 


Each term in the sum is the product of the value of the function at a given point and the length of an interval. 
Consequently, each term represents the (signed) area of a rectangle with height f (t;) and width х; + 1 — xj. 


4% 


Loosely speaking, the Riemann integral is the limit of the Riemann sums of a function as the partitions get 
finer. If the limit exists then the function is said to be integrable (or more specifically Riemann-integrable). 


“This part is sometimes referred to as the first fundamental theorem of calculus.” (Apostol 1967) 
Let f be a real-valued function defined on the interval [a, b]. Let F be the function defined by 
x 
F(x) = | fdt 
a 


The F is uniformly continuous on [a, b] and differentiable on the open interval (a,b) and 


F'(x) = f(x) 


For all x in (a, b) so F is an antiderivative of f. 


It’s Corollary, 


b 
| f'GOdx = F(b) — F(a) 


a 
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The two-dimensional case we have 


An n X m matrix, is a matrix with n rows and m columns A can be represented as 


йт? 


If A is an m x n matrix and B an n x p matrix, the matrix product С = AB is defined to be the m x p 
matrix such that 
n 

су = ag bi; + abo... фаућуј = У agb; for i = 1,...,m and j = 1,...,p. 

k=1 
That is, the entry c;; of the product is obtained by multiplying term-by-term the entries of the ith row of 
A and the jth column of B, and summing these n products. In other words cj; is the dot product of the 
ith row of A and the jth column of B. Thus the product of AB is defined if and only if the columns of 


A equals the number of rows in B. 
The identity matrix for n X n matrices is given by I, AI = А, I? = I, e.g. 159 = К | зат 


1 0 0 
0 1 0 
0 0 1 


Given an arbitrary matrix B = |: 2 
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Matrix Multiplication 
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The Euclidean group is the group of (Euclidean) isometries of a Euclidean space E”; that is, the 
transformations of space that preserve the Euclidean distance between any two points, that is the group of 
Euclidean transformations. The group depends only on the dimension n of space and is denoted [SO (n) for 
inhomogeneous special orthogonal group. The Euclidean group E (n) comprises all translations, rotations, and 
reflections of IE" ; and arbitrary finite combinations of them. The Euclidean group can be seen as the 
symmetry group of the space itself and contains the group of symmetries of any figure (subset) of that space. 
The continuous trajectories in /SO (3) play an important role in classical mechanics, because they describe the 
physically possible movements of a rigid body in three-dimensional space over time. One takes f (0) to be the 
identity transformation / of E? which describes the initial position of the body. The position and orientation of 
the body at any later time t will be described by the transformation f (t). Since f (0) = I is in E*(3), the same 
must be true of f (t) for any later time. For that reason, the direct Euclidean isometries аге also called "rigid 
motions". 


Euclidean metric transformations preserve invariant size of the figure, they only change the figure position in 
space. Relation between coordinates of the original and its image in the given transformation is expressed by 
equations of the transformation. The matrix of linear transformations is the matrix of the system of equations 


А = [х, у, 2,1] 
(х,у,2,1) > (х',у',2',1) 
х' = f(x,y,z) 

y = д(х,у,7) 


2' = h(z,y,z) 
А'=А.Т 


The identity is given by A’ = А · 1, Т = 1, = 


A'-A.T] 


For example, symmetry with respect to the coordinate plane л = xy, can be given by A’ = ATsr 


Rotation around a coordinate axis 2 by angle ф, A’ = А · То, 
cos@ sing 0 0 
T, =|7 sing соѕф 0 0 
E 0 0 10 
0 0 0 1 
х' -хсовф- уѕіпф, 
У => 


Translation by a vector у = (m, n, p, 0) 


= [х, у, 2,1] =[x+myt+n,z+p,1] 


1 
0 
0 


0 
x'=x+m 
y'=y+n 
z =z+p 


The general affine transformation is given by, |T,| = 0 


ТА = 


х' =ax+dy+gz+m 

y'=bx+ey+hz+n 

zi=cxt+fytiz+p 
h'=s 


А = (x,y,Zz,1) > А' (у 2 ВК) =A -Ta == 
(totem bx +ey + n cx+fy+iz+p 1) 
= s | :Ы —— ed 


, , 


5 5 5 


The coefficients a through to i, m n, p, s = 0 have а clear geometric interpretations: 


a,e,i are nonzero scales on coordinate axis x,y,z 

5 coefficient of scaling to the origin 

m,n, p — coordinates to the translation vector 

b, c, d, f, g, h — coefficients to the general affine transformation. 


Affine transformations do not preserve invariant the size of line segments and angles. 
Special subsets of affine transformations form similarities, preserving invariant the size of angles. 
Any Euclidean transformation is an affine transformation. 


The Galilean Transformation is given by: 


у 
"-2 


t'=t 
Note that the last equation holds for all Galilean transformations up to addition of a constant and expresses the 
assumption of a universal time independent of the relative motion of different observers. 


In the one-dimensional case we can show that newton’s second law is invariant under a Galilean 


transformation. 
Б d^x 
=т—— 
dt? 
The Galilean transformation in one dimension, with velocity in the positive x direction given by: 


x=x'+vt' 
t=t' 


y! 
Also note quickly rearranging this formula gives us the definition of velocity х = = у, that is, the change in 


distance over some global time that has taken place is our velocity. 


Let F' be the force in the primed coordinates. 


d’x d? RU d? x' 

de а” ар 

Hz d?x' Е а2х = 
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Note that m does depend on velocity relativistically, so the 2nd law and/or notion of invariance must be 
changed in relativity. In the relativistic case under the Lorentz transformation time is dependent on velocity 
and position and with reference to the speed of light. 


Again, in the one-dimensional case with velocity in the x direction. 


х' = y(x — vt) 
vx 
tS У (t = —) 
1 
where, y= 
172 
тыш 


When speed у is much smaller Шап с, the Lorentz factor is negligibly different from 1, but as у approaches с, 
ү grows without bound. The value of у must be smaller than с for the transformation to make sense. 


A Lorentz boost, which is simply a Lorentz transformation that doesn’t involve rotation, in direction n with 
magnitude v 


t' 2y(t—vn:x) 
х' =x + (у – 1)(Х · па — ут 


Assume two inertial reference frames (t, x, y, Z) and (t’, x’, y’,z’) and two points Р, and P}, the Lorentz group 


is the set of all transformations between the two reference frames that preserve the speed of light propagating 
between the two points. 


c* (At)* — (Ax')? — (ду)? — (42)? = c? (At)? — (Ах)? — (Ду)? — (Az)? 


In matrix form these are all the linear transformations A’nA = n 


ct’ -1 00 0 ct 
МИН | _|0 100 _ |x 
Х = |у| "7|o у у |у 
2' 0 001 7 


The spacetime interval takes the form 
X- X= XTX = X'TyX' 
And is invariant under the Lorentz transformation 


Х' = AX 
Where A is a square matrix which can depend on parameters. 


i. Light 


I'm not sure what light is. But I know I can see something over there. I know what a nice sunset looks like. I 
know what it feels like to stumble in darkness. Something is over there, and I am here. The relationship 
between those two things, and really that statement itself, is a definition of light, the photon, the emission and 
absorption of a quanta of energy between two particles. But how does the universe choose which particles to 
interact with which where and when. And what too. 


i. Gauss’s Law (for electricity) 


In its integral form, it states that the flux of the electric field out of an arbitrary closed 
surface is proportional to the electric charge enclosed by the surface, irrespective of 
how that charge is distributed. Even though the law alone is insufficient to determine 
the electric field across a surface enclosing any charge distribution, this may be 
possible in cases where symmetry mandates uniformity of the field. Where no such 
symmetry exists, Gauss's law can be used in its differential form, which states that the 
divergence of the electric field is proportional to the local density of charge. 


Gauss’s law may be expressed 
Ф Е = 9 
60 
Where Фу is the electric flux through а closed surface S, enclosing апу volume V, О 
is the total charge enclosed within V, and e, is the permittivity of free space. 


Фр = фЕ dA 
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Flux describes any effect that appears to pass or travel through a surface. 


ii.i) Flux 


First the flux at a single point is given 


J mom GM 
A 
Where j is the given flux, q is the physical quantity that flows, t is time and A 


is area, that is how much does the quantity flow over time over across an area. 


Flux as a scalar field defined a function of a set of point on a surface 


"ano 
TON Ар) 


dq 
1(A,p) = = Asp) 
As before, the surface is assumed to be flat, and the flow is assumed to be 
everywhere perpendicular to it. However, the flow need not be constant. q is 
now a function of p, a point on the surface, and A, an area. Rather than 
measure the total flow through the surface, q measures the flow through the 
disk with area A centered at p along the surface. 


Finally, flux as a vector field: 


Цр) 


і) = AQ) 


d 
I(A,p) = BiggestVector (5 = (A, p, 2 


In this case, there is no fixed surface we are measuring over. q is a function of 
a point, an area, and a direction (given by a unit vector й), and measures the 
flow through the disk of area A perpendicular to that unit vector. I is defined 
picking the unit vector that maximizes the flow around the point, because the 
true flow is maximized across the disk that is perpendicular to it. The unit 
vector thus uniquely maximizes the function when it points in the "true 
direction" of the flow. 


If the flux j passes through the area at an angle Ө to the area normal fi, the the 
dot product 
ј fi jcos0 


For vector flux the surface integral of j over a surface S, gives the proper 
flowing per unit through the surface 


dq 
тш |. AdA 
dt f; с 


ii.ii) Inductance 


Inductance is the tendency of an electrical conductor to oppose a change in 
the electric current flowing through it. The electric current produces a 
magnetic field around the conductor. The magnetic field strength depends on 
the magnitude of the electric current and follows any changes in the 
magnitude of the current. From Faraday's law of induction, any change in 
magnetic field through a circuit induces an electromotive force (EMF) 
(voltage) in the conductors, a process known as electromagnetic induction. 


Inductance is measured in Henry’s that have the SI base units 


kg -m? 
7252.42 
Where А is ampere, One ampere is equal to 1 coulomb (С) moving past a 


point per second and is defined in terms of the elementary charge e, at about 
6.241509 x 1018e The elementary charge, usually denoted by e, is a 


fundamental physical constant, defined as the electric charge carried by a 
single proton or, equivalently, the magnitude of the negative electric charge 
carried by a single electron, which has charge —1 e 


Two thin, straight, stationary, parallel wires, a distance r apart in free space, 
each carrying a current I, will exert a force on each other. Ampére's force law 
states that the magnetic force Fm per length L is given by 


(4: ) | 
2 Е Ho А 
а шины ы 

m L 2п т 


The vacuum magnetic permeability is given by, 


H 
Ho = 4T X 1077 — = 1.256...x 1076 kg. т.572. A7? 


1 
E0 =—> 
9 noc? 
By the divergence theorem, Gauss's law can also be written in the differential form 
p 
У-Еш- 
60 


Where У · Е is the divergence of the electric field, and р is the charge density, ог 
charge per unit volume. 


Илу) Divergence 


In cartesian coordinates the divergence of a continuously differentiable vector 
field F = Бі ФЕЈ + Ek 15 
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Vector fields are often illustrated using the example of the velocity field of a 
fluid, such as a gas or liquid. A moving liquid has a velocity—a speed and a 
direction—at each point, which can be represented by a vector, so that the 
velocity of the liquid at any moment forms a vector field. Consider an 
imaginary closed surface S inside a body of liquid, enclosing a volume of 
liquid. The flux of liquid out of the volume at any time is equal to the volume 
rate of fluid crossing this surface, i.e., the surface integral of the velocity over 
the surface. 


Since liquids are incompressible, the amount of liquid inside a closed volume 
is constant; if there are no sources or sinks inside the volume then the flux of 
liquid out of S is zero. If the liquid is moving, it may flow into the volume at 
some points on the surface S and out of the volume at other points, but the 
amounts flowing in and out at any moment are equal, so the net flux of liquid 
out of the volume is zero. 


However if a source of liquid is inside the closed surface, such as a pipe 
through which liquid is introduced, the additional liquid will exert pressure on 
the surrounding liquid, causing an outward flow in all directions. This will 
cause a net outward flow through the surface S. The flux outward through S 
equals the volume rate of flow of fluid into S from the pipe. Similarly, if there 
is a sink or drain inside S, such as a pipe which drains the liquid off, the 
external pressure of the liquid will cause a velocity throughout the liquid 


directed inward toward the location of the drain. The volume rate of flow of 
liquid inward through the surface S equals the rate of liquid removed by the 
sink. 


If there are multiple sources and sinks of liquid inside S, the flux through the 
surface can be calculated by adding up the volume rate of liquid added by the 
sources and subtracting the rate of liquid drained off by the sinks. The volume 
rate of flow of liquid through a source or sink (with the flow through a sink 
given a negative sign) is equal to the divergence of the velocity field at the 
pipe mouth, so adding up (integrating) the divergence of the liquid throughout 
the volume enclosed by S equals the volume rate of flux through S. This is the 
divergence theorem. 


The divergence theorem is employed in any conservation law which states 
that the total volume of all sinks and sources, that is the volume integral of the 
divergence, is equal to the net flow across the volume's boundary. 


Suppose V is a subset of IR” (in the case n = 3, V represents a volume in 
three-dimensional space) which is compact and has a smooth boundary S. 
If F is a continuously differentiable vector field defined on a neighbourhood 


of V then, 
|| (V. Е)ау = је й)а$ 


The integral form of Gauss's law is: 


фаал == 
5 бр 


for any closed surface 5 containing charge О. Ву the divergence theorem, this 


equation is equivalent to: 
| | | У · Еау = = 
60 
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for any volume V containing charge О. By the relation between charge апа 
charge density, this equation is equivalent to: 


ІШ e 


for any volume V. In order for this equation to be simultaneously true for 
every possible volume V, it is necessary (and sufficient) for the integrands to 
be equal everywhere. Therefore, this equation is equivalent to: 


Еее 
Eo 


The is the first of maxwells equations. 
Gauss’s Law for Magnetism or the absence of free magnetic poles. 
The differential form for Gauss's law for magnetism is: 


У-В-0 


iii. 


where B denotes the magnetic field. The integral form of Gauss’s law for magnetism 
states 


Фр = фв-45-о 
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The law in this form states that for each volume element in space, there are exactly the 
same number of "magnetic field lines" entering and exiting the volume. No total 
"magnetic charge" can build up in any point in space. For example, the south pole of 
the magnet is exactly as strong as the north pole, and free-floating south poles without 
accompanying north poles (magnetic monopoles) are not allowed. In contrast, this is 
not true for other fields such as electric fields or gravitational fields, where total 
electric charge or mass can build up in a volume of space. It is equivalent to the 
statement that magnetic monopoles do not exist. 


The magnetic field B can be depicted via field lines (also called flux lines) – that is, a 
set of curves whose direction corresponds to the direction of B, and whose areal 
density is proportional to the magnitude of B. Gauss's law for magnetism is equivalent 
to the statement that the field lines have neither a beginning nor an end: Each one 
either forms a closed loop, winds around forever without ever quite joining back up to 
itself exactly, or extends to infinity. 


Faraday's Law (of Inductance) 
The most widespread version of Faraday's law states: 


The electromotive force around a closed path is equal to the negative of the time rate 
of change of the magnetic flux enclosed by the path. 


For a loop of wire in a magnetic field, the magnetic flux Ф» is defined for any surface 
X whose boundary is the given loop. Since the wire loop may be moving, we write 
X(t) for the surface, The magnetic flux is the surface integral, 


Фр = $ B(t) -dA 
x(t) 


iii.i) The Maxwell-Faraday Equation 
The Maxwell—Faraday equation states that a time-varying magnetic field 
always accompanies a spatially varying (also possibly time-varying), non- 
conservative electric field, and vice versa. The Maxwell—Faraday equation is 


УХЕ = 28 
~ Qt 


A vector field whose curl is zero is called irrotational. 


In three dimensions the curl is given by 


t f UE 
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iv. 


The curl of a vector field denoted V х F, is an operator that maps k finitely 
differentiable smooth functions to k — 1 finitely differentiable smooth 
functions. Unlike the gradient and divergence, curl as formulated in vector 
calculus does not generalize simply to other dimensions; some generalizations 
are possible, but only in three dimensions is the geometrically defined curl of 
a vector field again a vector field. This deficiency is a direct consequence of 
the limitations of vector calculus; on the other hand, when expressed as an 
antisymmetric tensor field via the wedge operator of geometric calculus, the 
curl generalizes to all dimensions. 


Ampere’s circuital law law with Maxwell’s addition. 


In 1820 Danish physicist Hans Christian Ørsted discovered that an electric current 
creates a magnetic field around it, when he noticed that the needle of a compass next 
to a wire carrying current turned so that the needle was perpendicular to the wire. He 
investigated and discovered the rules which govern the field around a straight current- 
carrying wire: 


The magnetic field lines encircle the current-carrying wire. 

The magnetic field lines lie in a plane perpendicular to the wire. 

If the direction of the current is reversed, the direction of the magnetic field 
reverses. 

The strength of the field is directly proportional to the magnitude of the current. 
The strength of the field at any point is inversely proportional to the distance of 
the point from the wire. 


This sparked a great deal of research into the relation between electricity and 
magnetism. André-Marie Ampére investigated the magnetic force between two 
current-carrying wires, discovering Ampére's force law. In the 1850s Scottish 
mathematical physicist James Clerk Maxwell generalized these results and others into 
a single mathematical law. The original form of Maxwell's circuital law, which ће 
derived as early as 1855 in his paper "On Faraday's Lines of Force"[9] based on an 
analogy to hydrodynamics, relates magnetic fields to electric currents that produce 
them. It determines the magnetic field associated with a given current, or the current 
associated with a given magnetic field. 


The original circuital law only applies to a magnetostatic situation, to continuous 
steady currents flowing in a closed circuit. For systems with electric fields that change 
over time, the original law (as given in this section) must be modified to include a 
term known as Maxwell's correction. 


The integral form of the original circuital law is a line integral of the magnetic field 
around some closed curve C (arbitrary but must be closed). The curve C in turn 
bounds both a surface S which the electric current passes through (again arbitrary but 
not closed—since no three-dimensional volume is enclosed by S) and encloses the 
current. The mathematical statement of the law is a relation between the circulation of 
the magnetic field around some path (line integral) due to the current which passes 
through that enclosed path (surface integral). 


In terms of total current, (which is the sum of both free current and bound current) the 
line integral of the magnetic B-field (in teslas, T) around closed curve C is 
proportional to the total current lenc passing through a surface S (enclosed by C). In 
terms of free current, the line integral of the magnetic H-field (in amperes per metre, 
A · тг) around closed curve C equals the free current If enc through a surface S. 


Using B-field and total current 
The Integral form is given, 


$ B - d1 = po PJ - dS = polene 
C S 


The differential form is given 
Vx B= ш] 


J is the total current density (in amperes per square meter. А · m 2) 
$- is the closed line integral around closed curve C 


ff, denotes a surface integral over the surface S bounded by the curve C 


dl is an infinitesimal element (a differential) of the curve C (i.e. a vector with 
magnitude equal to the length of the infinitesimal line element, and direction given by 
the tangent to the curve C) 

dS is the vector area of an infinitesimal element of surface S (that is, a vector with 
magnitude equal to the area of the infinitesimal surface element, and direction normal 
to surface S. The direction of the normal must correspond with the orientation of C by 
the right-hand rule), see below for further explanation of the curve C and surface S. 


There are two important issues regarding the circuital law that require closer scrutiny. 
First, there is an issue regarding the continuity equation for electrical charge. In vector 
calculus, the identity for the divergence of a curl states that the divergence of the curl 
of a vector field must always be zero. Hence, 


V-(VxB)=0 
And so Ampere’s circuital law implies, 
V-J=0 


i.e. that current density is solenoidal (solenoidal by definition means V v = 0.) 
But in general, to whatever it means that reality follows the continuity equation for 
electric charge is given, 


Which is nonzero for a time varying charge density. 
iv.i) Continuity equations 


A continuity equation or transport equation is an equation that describes the 
transport of some quantity. It is particularly simple and powerful when 
applied to a conserved quantity, but it can be generalized to apply to any 
extensive quantity. Since mass, energy, momentum, electric charge and other 
natural quantities are conserved under their respective appropriate conditions, 
a variety of physical phenomena may be described using continuity equations. 


iv.i.i) With regards to flux 


A continuity equation is useful when a flux can be defined. 


Фр = E -a4 
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And more generally, 
Flux (the rate q flowing through the imaginary surface S) 


- [| а 


e.g. electric flux 


Figure I Illustration of how the fluxes, or flux densities, jq and Jz of a 
quantity q pass through open surfaces S4 апа S5. (vectors Sqand S represent vector 
areas that can be differentiated into infinitesimal area elements). 


The integral form of the continuity equation states that: 


The amount of q in a region increases when additional q flows inward 
through the surface of the region, and decreases when it flows 
outward; The amount of q in a region increases when new q is created 
inside the region, and decreases when q is destroyed; 


That is, 


+ bj dS = У 
atf 


У is the net rate q is being generated inside a volume V per unit time when q 
is being it is called a source of q, and it make X positive, when q is being 
destroyed it is called a sink of q, and it makes X more negative. This term is 
sometimes written dq/dt| gen or the total charge of q from it generation or 
destruction inside the control volume. 
By the divergence theorem, a general continuity equation can also be written 
in differential form, 

др 
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Regarding the issue the propagation of electromagnetic waves, in free space where, 


Ј= 0 


VxB=0 
i.e. that the magnetic field is irrotational but to maintain consistency with the 
continuity equation for electric charge we must have, 


The circuital law implies 


To treat these situations, the contribution of displacement current must be added to the 
current term in the circuital law. 


James Clerk Maxwell conceived of displacement current as a polarization current in 
the dielectric vortex sea, which he used to model the magnetic field hydrodynamically 
and mechanically.[17] He added this displacement current to Ampére's circuital law at 
equation 112 in his 1861 paper "On Physical Lines of Force”. 


Displacement Current 


In free space, the displacement current is related to the time rate of change of 
electric field. 


The Model 


Given a simple harmonic oscillator such as a spring: 


time 
me, 


The work done by a conservative force, such in the case of the spring, the gravitational force, is 
W = —AU 
where AU is the change in the potential energy associated with the force. 
Kinetic Energy + Potential Energy = Total Energy 
K.E+P.E=E 


The change in kinetic energy is equal to the work done associated with the force and displacement by newtons second law 
of motion. 


AK =W 


dv dx 
AK = | F dx = | madx =m | Z ах =m | ағ-т| vav 
dt dt 


So we consider the change in kinetic energy, to be the force exerted at each point along some displacement, that is 
proportional to its mass and acceleration, proportional to its momentum and its velocity. 


v 1 1 
AK — m | vdv = —mv? — – тд 
" 2 2 


Naturally, the kinetic energy of an object at rest should be zero. Thus, an object's kinetic energy is defined mathematically 
by the following equation... 


Classically Kinetic Energy 


Classically Potential Energy follows in the same fashion 
Р.Е = | Е-ах 
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P.E = =kx? 
2 
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=kx? = E 
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Momentum in quantum mechanics is given by 
ШЕ 
Pia 
Swapping this into our equation gives us 
h^ 0^ 1 kx? 
2móx 2^ 
We will have any arbitrary potential U(x) 
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We introduce the wave function w(x, t) 
The Hamiltonian form of the Schrodinger equation is given: 
Hy (x,t) = Еф(х, 0) 


That is, 


2m dx? 


MD 
+ P.E(x,t) |w(x,t) = Ew(x,t) 


(Kinetic Energy + Potential Energy)ParticleState = Energy x ParticleState 
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Which becomes. 
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Which is, 
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Thus, 
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Say we have a 2x2 space. At each point in space there is a complex number. There are 2? states. 
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If A is n X n, B is m x mand Ig denotes the k x k identity matrix, then we can define a “Kronecker sum” © 
АФВ = AGL, + 1,98 


The Kronecker product of matrices corresponds to the abstract tensor product of linear maps. Specifically, if the vector 
spaces V, W, X and Y have bases (vi, ..., Vm}, (Ил, ..., 24) (21, = Xa} and (ул, ..., Ve}, respectively, and if the matrices A 
and B present linear transformations S : У — X and T : W —~Y, respectively in the appropriate bases, then the matrix 
A GB represents the tensor product of two maps, 5 ӨТ: У ФУ — ХОУ with respect to the basis 

(и, 8w, v4 и, ... VW, ..., v, и} of V@W and the similarly defined basis of X&Y with the property that 

A GB(v;&w;) = (Av,)@(Bw,), where i and j are integers in the proper range. 
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It is going to be our attempt to solve the linear equation: 
Аф(х, =i+1) = By(x,t =i) 
First, we will introduce a complex valued wave function. That is 


s € W(x, t), sec 
Our space is some finite subset of Z? x = (x,y,z) | x,y,z € 73. Our space has a length L in each direction, for 
simplicity we consider a square space where = <x< 5, There are N X N X N points. 


Where Every point is a complex number. We recall the isomorphism of the structure of C to IR? 


We have two real spaces that are independent. 


We will start with a simple gaussian distribution for our particle state у(х, = 0) 


Vx, 0) = e2? 


We can modify the standard deviation o and centre и of the gaussian distribution 


Position Space 
7 Momentum Space 


2d-Gaussian distribution 
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We add a complex part, we introduce phase velocity into the gaussian, this part exists in the momentum space. 


w(x, 0) = er . g iPox 
For w(x,t = i + 1), from an initial wave function w(x, t = 0), 
Given, 
H=K.E+P.E 
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There is ап М x М tridiagonal matrix T, The matrix T represents the discrete Laplacian operator, a central 
difference approximation to the second derivative. 
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This matrix is then scaled by -— 


ae ae which represents the kinetic energy term in the finite difference 
approximation. 


We use Kronecker products to extend the 1D kinetic operator T to 3D space. 
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For every i, j and k position in our space we have some potential, 
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For every i, jand Кіп N there is an element of the potential U(ij.k) and w such ару =i-N-N+j:-N+k 
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We have the diagonal matrix 


Up 0 0 
10 U, : 
PES 1 0 
0 0 UW 
Where each diagonal entry is an element of our 3D potential. 
Uu 0 => 0 
0 U, - : 
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We аге solving for x 
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We use an iterative method to solve such as SciPy's GCROT 
w(t + 1) = x = GCROT(A,b) 
We will first explore an algebraically similar method. The conjugate residual method. 


The conjugate gradient method is an algorithm for the numerical solution of systems of linear equations. 


Suppose we want to solve the systems of linear equations 
Ax = b 


For the vector x. where the known n x n matrix А is symmetric i.e АТ = А, positive-definite (i.e x7 Ax > 0 for all 
nonzero vector x € IR") and real, and b is known as well. We denote the unique solution of this system x,. 


We say that two non-zero vectors u and v are conjugate (with respect to A) if 
ТА) = 
u Av = 0 


Since A is symmetric and positive definite, the left-hand side defines an inner product, 
ul Av = (u,v), = (Au, v) = (и, ATv) = (и, Av) 


An inner product is a generalization of the dot product. In a vector space, it is a way to multiply vectors together, with the 
result of this multiplication being a scalar. 


A function (.,.) : R” x R” > R is an inner product on R” if an only if there exists a symmetric positive-definite matrix 
M such that (x, y) = x1 My for all x, y € R”. If M is the identity matrix, then (x, y) = х My is the dot product. 


For R?, where M is positive-definite (which happens if an only if det M = ad — b? > 0 and one/both diagonal elements 
are positive) then for any x = [xi, x]7, y = [yi y2]" € R?, 


a b 


У. 
b d b.] = axyy, + Dx y2 + bx5y, + аху 


(x,y) = x! My = [xa] | 
Every inner product space on R? is of this form (where ђе Ка > 0 апаа > 0 satisfy ad > 52) 
The general form of an inner product on C” is known as the Hermitian form and is given by: 
(x,y) = ytMx = x'My 
The inner product of a pair of complex numbers z and w is the product of z with the complex conjugate of w: 
(z,w) = zw 


A second example is the space C? whose elements are pairs of complex numbers 2 = (21,2). Then the inner product of z 
with another such vector w = (w4, w3) is given by, 


(z, w) = ZW + Z2W; 


A Hilbert space is a real or complex inner product space that is also a complete metric space with respect to the distance 
function induced by the inner product. To say that a complex vector space H is a complex inner product space means that 
there is an inner product (x, y) associating a complex number to each pair of elements x, y that satisfies the following 
properties, 


1. The inner product is conjugate symmetric; that is, the inner product of a pair of elements is equal to the complex 
conjugate of the inner product of the swapped elements. 


(у,х) = (х,у), 
Importantly this implies that (x, x) is a real number. 


2. The inner product is linear in its first argument. For all complex numbers a and b, 
(ax, + рх, y) = a(x, y) + b(x2,y) 


3. The inner product of an element with itself is positive-definite: 
(х,х) > 0 ifx#0 
(х,х) = 0 ifx=0 
It follows from 1 and 2 that a complex inner product is antilinear (conjugate linear) in its second argument, meaning that, 


(x. ay, + by;) = G(x, y,) + b(x, y2) 


Suppose V is a vector-space. A function f that converts ordered pairs of vectors into scalars is a bilinear form if it satisfies 
the rules: 


e f(utv,w) = (шм) + f(v,w) 
• f(kuv) = k - fv) 
• fwv+w)=fuv)+fuw) 
e f(wkv) = К: flv) 


If f is a bilinear form, it is symmetric if it satisfies: 


• fv) = (о, и) 


If f is a symmetric bilinear form, it is positive definite if it satisfies: 
e 6f(u,u) > 0,where u = 0 


An inner product is a positive-definite symmetric bilinear form (While a real inner product (x, y) is bilinear (linear in 
each variable separately), there is a slight difference for complex inner products, A complex inner product (x|y) is linear 
in x and is conjugate linear in y). 


An inner-product space is a vector space with an inner product. By positive-definiteness, the norm |v| of a vector v in an 
inner-product space can be defined to be the non-negative square root of (v, v), ||v]| = y (v, v) 


The real numbers R are a vector space over R that becomes an inner product space with arithmetic multiplication as its 
inner product, 


(x.y):2 xy forx,y ER 
The complex numbers С are a vector space over C that becomes an inner product space with the inner product, 
(x,y):xy forx,y EC 


More generally the real n-space R” with the dot product is an inner product space, an example of a Euclidean vector 


space. 
Хі 
ІН 


Two vectors u and v are conjugate if an only if they are orthogonal with respect to the inner product 


У 


" 


n 
| =хТ у= Xiyi 
Yn i=1 


ul Av = (u,v), = (Au, v) = (и, ATv) = (и, Av) 
Being conjugate is a symmetric relation: if и is conjugate to у, then у is conjugate to u. 


Suppose that 


Р = Ру Pr} 


Is a set of n mutually conjugate vectors with respect to A, i.e. pi Ap; = 0 for alli = j. Then P forms a basis for IR", and 
we may express the solution x, of Ax — b in this basis: 


n n 


х, = У ар; > Ах, = > aj Api 


i=1 i=1 


Left multiplying the problem Ax = b with the vector pz yields, 


n n 
pkb = py Ax, = > ар Ap; = 2 ai(Pr Pila = ак (Pro Pila 
i=1 i=1 
And so, 
_ Феба 
pie cle 


7 (Pr Реја 


This gives the following method for solving the equation Ax = b: find a sequence of n conjugate directions, and then 
compute the coefficients ак. 


As an iterative method, if we choose the conjugate vectors p; carefully, then we may not need all of them to obtain a 
good approximation to the solution x,. So, we want to regard the conjugate gradient method as an iterative method. This 
also allows us to approximately solve systems where n is so large that the direct method would take too much time. 


We denote the initial guess for x, by xo, we can assume without the loss of generality that xy = 0. Starting with xy we 


search for the solution and in each iteration, we need a metric to tell us whether we are closer to the solution x,. This 
comes from the fact x, is also the unique minimizer of the following quadratic function 


1 
ја) = УХ Ax — xT b, x ЄЙ" 


The existence of a unique minimizer is apparent as its Hessian matrix of second derivatives is symmetric positive-definite 
H(f@))=A4 
And that the minimizer Vf (x) solves the initial problem follows from its first derivative, 
Vf (x) = Ax Б 


This suggests taking the first basis vector pọ to be the negative of the gradient of f at x = хо. The gradient of f equals 
Ax — b. Starting with an initial guess Хо, this means we take pọ = b — Ах. The other vectors in the basis will be 
conjugate to the gradient, hence the name conjugate gradient method. Note that pọ is also the residual provided by this 
initial step of the algorithm. 


Let ту be the residual at the kth step 

Т, = Ь = Ах, 
As observed above, ту, is the negative gradient of f at хк, so the gradient descent method would require to move in the 
direction of rg. Hence, however, we insist that the directions pg must conjugate to each other. A practical way to enforce 
this is requiring that the next search direction be built out of the current residual and all previous search directions. This 


conjugate constraint is an orthonormal-type constraint and hence the algorithm can be viewed as an example of a Gram- 
Schmidt orthonormalization. This gives the following expression: 


ick 
Following this direction, the next optimal location is given by, 


Xk+1 = Хк + акрк 
With 


рк (b — Axx) рк Атк 
ак = T E 
Di Арк Di Арк 


Where the last equality follows from the definition of rg. The expression for a, can be derived if one substitutes the 
expression for x444 into f and minimizing it with respect to ак 


f rar) = / (хк + акрк) = д(ак) 


Т. 
, ! pix (b — Axx) 
g' (ak) = 02 
i Ы p} Apr 


The above algorithm gives the most straightforward explanation of the conjugate gradient method. Seemingly, the 
algorithm as stated requires storage of all previous searching directions and residue vectors, as well as many matrix— 
vector multiplications, and thus can be computationally expensive. However, a closer analysis of the algorithm shows 
that r; is orthogonal to 7; , i.e. "m = 0, for i + j. And p; is A-orthogonal to p; , i.e. pi Apj = 0, for i + j This can be 
regarded that as the algorithm progresses, p; and r; span the same Krylov subspace. Where 7; form the orthogonal basis 


with respect to the standard inner product, and р; form the orthogonal basis with respect to the inner product induced 
by A. Therefore х, can be regarded as the projection of x on the Krylov subspace. 


That is if the conjugate gradient method starts with Хо = 0, then 
хк = argminycan{(x — y)! A(x — y) : y Є span{b, Ab, ... АЗЫ) 


The algorithm below is detailed for solving Ax — b where A is real, symmetric, positive definite matrix. The input vector 
хо can be an approximate initial solution or 0. It is a different formulation of the exact procedure described above. 


то:= b — Axo 
If то is sufficiently small, then return Хо as the result, else 


Do ‘= To 
k :=0 
repeat 
ттк 
ар '= -T 
PkAPk 


Хк+1 = Xk + акрк 
Tk+1 1— Tk — AAP, 
If y 44 is sufficiently small, then exit loop, else 


T 
B .— Патка 
к = HUS 


"тұ 
Dii "ұза + Ёкрк 
к:= К+1 


end repeat 
return Хк+1 as the result 


The conjugate gradient method with a trivial modification is extendable to solving, given complex-valued matrix A and 
vector b, the system of linear equations Ax — b for the complex-valued vector x, where A is Hermitian (i.e., A' 2 A) and 
positive-definite matrix, and the symbol ' denotes the conjugate transpose. The trivial modification is simply substituting 
the conjugate transpose for the real transpose everywhere. 


Appendix. 


ii. Series 


The exponential function is given as 


a x" x? x3 
ех = У Esite а + 
And sine as = 
| ee (—1)"х2"+1 x3 x? x? 
00 = A улу C 378 7 


These are all power series and are also all examples of Taylor Series. 


A Fourier series is a continuous, periodic function created by a summation of harmonically related sinusoidal 
functions. 


uae) | (Аһ соз (2х) + В, sin ( (277 x) јео Cne my 


n=1 n--N 


The first four partial sums of the Fourier series for a square wave. As more harmonics are added, the partial 
sums converge to (become more and more like) the square wave: 


| М / 1 


The Bloch or Poincare sphere representation for SU(2) where the base manifold is the S? sphere while the fibre 
is given by the U(1) phase at each point on that sphere. Together we have the finer bundle SU(2) = S? x U(1). 
Thus, S? the 3-sphere, is a principal U(1)-bundle over the 2-sphere 


The Bloch sphere is a geometrical representation of the pure state space |y} of a two-level quantum 
mechanical system. The Bloch sphere is a unit 2-sphere, with antipodal points corresponding to a pair of 
mutually orthogonal state vectors. The north and south poles of the Bloch sphere are typically chosen to 
correspond to the standard basis vectors |0) and |1), which in turn might correspond e.g. to the spin-up and 


spin-down states of an electron. The choice however is arbitrary. Given an orthonormal basis, any pure state 
|9) of a two-level quantum system can be written as a superposition of the basis vectors |0) and |1), where the 


|0) 


|0) + 11) 


10) +11) 
v2 


Ө ЖРК. 
ly) = cos = |0) + e*sin = |1) 


coefficient of (or contribution from) each of the two basis vectors is a complex number. This means that the 
state is described by four real numbers. 


(pip) = 1, or equivalently || |W) 12= 1, Given this constraint, that the pure state is unitary, we can write |p) 
using the following representation, 


фу = cos(0/2)|0) + e'?sin(8/2)|1) = соз (0 /2)|0) + (cos + isin$)sin(0/2)|1),where0 < 0 < пап40 < ф < 2п 


The representation is always unique, because even though the value of ф is not unique when |) is one of the 
states |0) ог |1), the point represented by бапаф is unique. 


i. Light at a point in Space as a complex number with real electric component and 
imaginary magnetic component. (The Silberstein Vector) 


Heinrich Martin Weber published the fourth edition of "the partial differential equations of mathematical 
physics according to Riemann’s lectures” in two volumes in 1900 and 1901. Weber demonstrated how to 
consolidate Maxwell’s equations without charges or currents, 


i O(E + iM) 


curl(E + iM) = 5 aE 


This would be independently rediscovered in 1907 by Ludwik Silberstein. 


Given an electric field E and a magnetic field B defined on a common region of spacetime, the Riemann- 
Silberstein vector is 


F = Е +ісВ 


Where c is the speed of light (с = 3 · 108ms~*). It is isomorphic to the electromagnetic tensor F, a 2-vector 
used in the covariant formulation of classical electromagnetism. The Riemann-Silberstein vector is used as a 
point of reference in the geometric algebra formulation of electromagnetism. 


In a 1996 contribution to quantum electrodynamics Iwo Bialynicki-Birula used the Riemann-Silberstein vector 
as the basis for an approach to the photon, noting that it is a “complex vector-function of space coordinates 
and time that adequately describes the quantum state of a single photon.” To put the Riemann-Silberstein 
vector in contemporary parlance, a transition is made, “With the advent of spinor calculus that superseded the 


quaternionic calculus, the transformation properties of the Riemann-Silberstein vector have become even more 
transparent ... a symmetric second-rank spinor.” 


Maxwell’s four equations in vector calculus reduce to one equation in the algebra of physical space. 


(2 ЊЕ = 2(р 27) 


The Riemann-Silberstein vector is a complex unit vector, F € С, it is thus an element of the Clifford algebra, 
an algebra generated by a vector space with a quadratic form, that is a polynomial with terms all of degree 2 
and is. a unitary associative algebra, Cl, , (IR),a two dimensional algebra generated by a single basis vector that 
squares to -1, and is algebra-isomorphic to C (С = Cl,(C}). In contrast to the electron wave function the 
modulus square of the wave function of the photon (Riemann-Silberstein vector) is not dimensionless and 
must be multiplied by the “local photon wavelength” with the proper power to give dimensionless expression 
to normalise. 


ii. The sphere 


52, the 2-sphere, commonly the sphere, is the boundary of a 3-ball in 3-dimensional space. 
A point on the two sphere can be parameterised with 3 cartesian coordinates J x? + y? + 22 = 1 
Or by two angles, 0, ф, an angle in the x-y plane and one between the x-z, the angles of longitude and latitude. 


a = (віпбсое<ф, ѕіпӨѕіпф, cos0) = (х,у,2) 


It is known that among the n-spheres, only 59, Stand S? can admit a lie group structure, which is a highly non- 
trivial theorem. 


SU(2) = (5 Б) :a,b € C, |a|? + |b|? = 1} Is diffeomorphic to the 3-sphere 53, if we consider a, В as a pair 


in C? where a = a + bi and f = c + di then the equation |а]? + |b|? = 1 becomes a? + b? + c? + d? = 1. SU(2) 
is isomorphic group of versors. Quaternions of norm 1 are called versors since they generate the rotation group 
SO(3) 


The SU(2) matrix 


е p) 


-сійі a- bi GPC ER) 


can be mapped to the quaternion a1 + bi + cj + dk. This map is in fact the group isomorphism that relates 
SU(2) to versors. 


One basis for SU(2) is given by, t4 = aM al ‚5 = | 2: ‚з= an | These are related to the Pauli 


| 1 
matrices by, t; ө 201 
1 


The unitary group U(n) is а real Lie group of dimension n2. Тһе Lie algebra of U(n) consists of n x n skew- 
Hermitian matrices. 


In linear algebra it is often important to know which vectors have their directions unchanged by a linear 
transformation. An eigenvector or characteristic vector is such a vector. Thus an eigenvector v of a linear 
transformation T is scaled by a constant factor 2 when the linear transformation is applied to it Tv = Av. 


Geometrically vectors are multi-dimensional quantities with magnitude and direction, often pictured as 
"arrows". A linear transformation rotates, stretches, or shears the vectors it acts upon. Its eigenvectors аге 
those vectors that are only stretched with no rotation or shear. The corresponding eigenvalue is the factor by 
which an eigenvector is stretched or squished. 


Hermitian matrices are fundamental to quantum mechanics because they describe operators with necessarily 


real eigenvalues. An eigenvalue a of an operator A on some quantum state |Y) is one of the possible 
measurement outcomes of the operator, which requires the operators to have real eigenvalues. 


A Hermitian matrix or self-adjoint complex square matrix that is equal to its own conjugate transpose— that 
is, the element in the ith row and and jth column is equal to the complex conjugate of the element in the jth 
row and ith column for all indices i and j. 

e а; 


— ; Ke 
аң or in matrix form Apermitian €? A = А 


A Hermitian j— 


The group U(2) is the set of all 2 x 2 matrices such that U*U = I. 
Every matrix of the lie algebra has the form 
al + ро + 5502 + b303 
Where / is the 2 x 2 identity matrix, 61, 02, оз are the three Pauli matrices, and a4, b4, bz, b3 € R. 
This can be rewritten in the form 
al t b-.o 
Where b = (b4, bz, b4) € R? and o = (04, 02,03) 
U (2) = {exp(ial + ib - o)|a € Rand b є ЕЗ) 


Pauli matrices are Hermitian, involutory and unitary. Each Pauli matrix is Hermitian and together with the 
identity matrix 1 form a basis of the real vector space of 2 x 2 Hermitian matrices. This means that any 2 x 2 
Hermitian matrix can be written in a unique way as a linear combination of Pauli matrices, with all coefficients 
being real numbers. Pauli matrices span the space of observables of the complex 2 dimensional Hilbert space. 
The Pauli matrices after multiplication by i to make them anti-hermitian also generate transformation in the 
sense of lie algebras: the matrices io,, io;, ѓоз for a basis for the real lie algebra SU(2), which exponentiates to 
the special unitary group SU(2). The algebra generated by the three matrices c1, 05,05 are isomorphic to the 
Clifford algebra of R?, and the unitary associative algebra generated by 11,10, io; are isomorphic to that of 
quaternions. 

All three of the Pauli matrices can be compacted into a single expression: 


0; = у 


бук 18 the kronecker delta which equals +1 if j = К and 0 otherwise. 


The matrices are involutory: 

02 = д2 = оз = —i0 10,03 = 1; 
With the inclusion of the identity matrix, the Pauli matrices form an orthogonal basis of the Hilbert space of 
2 x 2 Hermitian matrices, J£; over R, and the Hilbert space of all complex 2 x 2 matrices over C. 


W. 
W Тапа W- bosons 1“ 
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2- bosonand y { B 

Mathematically, electromagnetism is unified with the weak interactions as a Yang-Mills field with an SU(2) x 
U(1) gauge group, , which describes the formal operations that can be applied to the electroweak gauge fields 


without changing the dynamics of the system. These fields are the weak isospin fields W1, W2, and W3, and 
the weak hypercharge field B. This invariance is known as electroweak symmetry. 


The generators of SU(2) and U(1) are given the name weak isospin T and weak hypercharge Y respectively. 
These then give rise to the gauge bosons that mediate the electroweak interactions — the three W bosons of 
weak isospin (Ил, W2, and W3), and the B boson of weak hypercharge, respectively, all of which are 
"initially" massless. These are not physical fields yet, before spontaneous symmetry breaking and the 
associated Higgs mechanism. 


In the Standard Model, the observed physical particles, ће W+ and 2, bosons, and the photon, are produced 
through the spontaneous symmetry breaking of the electroweak symmetry SU(2) x U(1), toU()em, 
effected by the Higgs mechanism (see also Higgs boson), an elaborate quantum-field-theoretic phenomenon 
that "spontaneously" alters the realization of the symmetry and rearranges degrees of freedom. 


The electric charge arises as the particular linear combination (nontrivial) of YW (weak hypercharge) and the 
T3 component of weak isospin (Q = T3 + 5 Yy)that does not couple to the Higgs boson. That is to say: the 


Higgs and the electromagnetic field have no effect on each other, at the level of the fundamental forces ("tree 
level"), while any other combination of the hypercharge and the weak isospin must interact with the Higgs. 
This causes an apparent separation between the weak force, which interacts with the Higgs, and 
electromagnetism, which does not. Mathematically, the electric charge is a specific combination of the 
hypercharge and T3. 


U(1)em (the symmetry group of electromagnetism only) is defined to be the group generated by this special 
linear combination, and the symmetry described by the U(1) em group is unbroken, since it does not directly 
interact with the Higgs. 


The above spontaneous symmetry breaking makes the W3 and B bosons coalesce into two different physical 
bosons with different masses — the 2 boson, and the photon (у), 


s xu a TE 
29 \=sinOy cos Ow) Из 


where Ow is the weak mixing angle. The axes representing the particles have essentially just been rotated, in 
е (W3, B) plane, by the angle Ow. This also introduces a mismatch between the mass of the 40 and the mass 
of the WE particles (denoted as mz and my, , respectively), 


Ту 


mz = —— 
2” cos Oy 


Тһе W; and W, bosons, in turn, combine to produce the charged massive bosons W+: 


wt=—w + iW) 
VZ 1 2 


Because the value of the mixing angle is currently determined empirically, in the absence of any superseding 
theoretical derivation it is mathematically defined as: 


My 
cos Q,, = — 
Mz 


The gamma matrices {y1, y?, уз, y^) are a set of 4 x 4 complex matrices which satisfy the defining anti- 
commutation relations that ensure they generate a matrix representation of the Clifford algebra Cl, з (К). When 
interpreted as the matrices of the action of a set of orthogonal basis vectors for contract variant vectors in 
Minkowski space, the column vectors on which the matrices act become a space of spinors, on which the 
Clifford algebra of spacetime acts, this in turn makes it possible to represent infinitesimal spatial rotations and 
Lorentz boosts. 


ү° is the time-like hermitian matrix. The other three are space-like, anti-Hermitian matrices. More compactly 
y? = о? @ L, andy! = io? ® сі, where & denotes the Kronecker product (Block matrix multiplication maybe 
more detail needed). 


The group SU(2) is the Lie group of unitary 2 x 2 matrices with unit determinant; its Lie algebra is the set of 
all 2 x 2 anti-Hermitian matrices with trace 0. The Lie algebra SU(2) is the 3-dimensional real algebra spanned 
by the set io, 

SU(2) = spans, 2,3) 

The lie algebra SU(2) is isomorphic to the lie algebra SO(3), which corresponds to the Lie group SO(3), the 
group of rotations in three-dimensional space. In other words, one can say that the io; are a realisation and in 
fact, the lowest-dimensional realisation) of infinitesimal rotations in three-dimensional space. However even 
though SU(2) and SO(3) are isomorphic as Lie algebras, SU(2) and SO(3) are not isomorphic as Lie groups. 
SU(2) is actually a double cover of SO(3), meaning that there is a two-to-one group homomorphism from 
SU(2) to SO(3). 


The real linear span of (1, 01,10. 104) is isomorphic to the real algebra of quaternions, H, represented by the 
span of the basis vectors {1, i, j,k}. The isomorphism from H to this set is given by the following map 


1» liv —0504 --ібі,) > —030, = —io,,k > —010; = —103 
The spin group is used to describe the symmetries of (electrically neutral, unchanged) fermions. It's 
complexification, Spin, is used to describe electrically charged fermions, most notably the electron. Strictly 


speaking, the spin group describes a fermion in a zero-dimensional space; however, space is not zero- 
dimensional, and so the spin group is used to define spin structures on pseudo-Riemannian manifolds. 


Each of the Hermitian Pauli matrices of spin— 5 particles has two eigenvalues, +1 and -1. The corresponding 


normalised eigen vectors are, 


1 
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a+ =|Ь%}, = A 
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By the postulates of quantum mechanics, an experiment designed to measure the electron spin on the x,y, or z 
axis can only yield an eigenvalue of the corresponding spin operator (S,,5y,orS,) on that axis, i.e. t or — - The 


quantum state of a particle with respect to spin can be represented by a two-component spinor: 


UM 


... lots to talk about spinor 
The Dirac equation is written in terms of a Dirac spinor field y taking values in a complex vector space 
described concretely as С“, defined on a flat space time R+’. 


In terms of a field y: 813 > C? , the Dirac equation is then, 
(ту"д, – mc)y(x) = 0 
(уу? = 21441, 
These matrices are explicitly realised under a choice of representation, Pauli matrices аге often chosen. 
ics | 0 z] 
ps -в 0 


Square-integrable function L? 


f : R > C square integrable = Í |f x) Рах < со 
Or more generally, as defined by the inner product, 
б.а) = | 700 водах 
А 

уућеге, 

f and g are square integrable functions, 

f (x) is the complex conjugate of f(x), 

since |a|? = a - а, square integrability is the same as (f, f) < оо 
An equivalent definition is to say that the square of the function itself (rather than of its absolute value) is 
Lebesgue integrable. For this to be true, the integrals of the positive and negative portions of the real part must 
both be finite, as well as those for the imaginary part. 
The space of square integrable functions is the L? space in which p = 2. 
The vector space IR", together with the dot product, is an example of an inner product space. In particular, IR 
itself (as a vector space), with multiplication, is an inner-product space. 


